Bipartite quantum states and random complex networks 
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We introduce a mapping between graphs and pure quantum bipartite states and show that the 
associated entanglement entropy conveys non-trivial information about the structure of the graph. 
Our primary goal is to investigate the family of random graphs known as complex networks. In 
the case of classical random graphs we derive an analytic expression for the averaged entanglement 
entropy S while for general complex networks we rely on numerics. For large number of nodes n 
we find a scaling S ~ clogra + g e where both the prefactor c and the sub-leading 0(1) term g e are 
a characteristic of the different classes of complex networks. In particular, g e encodes topological 
features of the graphs and is named network topological entropy. Our results suggest that quantum 
entanglement may provide a powerful tool in the analysis of large complex networks with non-trivial 
topological properties. 



Introduction. Complex networks are models of graphs 
that appear able to capture the phenomenology of a 
plethora of systems, from biology to the World Wide 
Web Q. The departure from regular lattices, the most 
common background geometry in solid state physics, al- 
lows for the rich static and dynamic behavior of com- 
plex networks. This is due to the simultaneous presence 
of a global compact structure and a sophisticated archi- 
tecture of interactions. By compact structures we refer 
to the typical small distance (with respect to a regular 
lattice) between nodes in the network. The complexity 
in the network architecture is manifest in the entangled 
pattern of links and paths that such objects display, see 
Figffl This architecture encodes a type of strong disorder 
that requires for its analysis some of the techniques devel- 
oped in Statistical Mechanics Q. There exist a plethora 
of phenomenological quantities that provide information 
on the architecture of a network: degree distribution, 
clustering coefficient, community structure measures and 
many others [![. In particular a few classical entropic 
measures have been introduced to describe the structure 
of complex networks 

In this letter we address the problem of the entropic 
analysis and discrimination of networks using quantum 
information tools, notably entanglement entropy. Entan- 
glement, a purely quantum measure of correlation, is one 
of the fundamental concepts in quantum information 0] ■ 
We provide a recipe, in a way the simplest possible one, 
to construct a pure bipartite quantum state for any given 
graph. This allows us to study entanglement properties 
of quantum states that are related to the topological fea- 
tures of the original graphs, and that are able to dis- 
tinguish between different complex network topologies. 
Although at the first sight it may seem a bit artificial to 
look for a graph-entropy measure in a quantum context 
the synergy between quantum information and complex 
network tools is not new. For example, in and @ the 
authors have discussed different interesting ways to asso- 





FIG. 1: (Color online) Random graphs belonging to different 
ensemble: Erdos-Renyi (top left corner), preferential attach- 
ment (top right corner), small- world (bottom left corner), and 
the same graph but with a different embedding (bottom right 
corner) . 



ciate graphs to quantum states and investigated in which 
sense complex networks may play a role in the quantum 
domain. All these constructions are then similar in spirit 
but substantially different from the present approach. 

The paper is organized as follows. We first describe 
the construction of the quantum bipartite network states. 
Then we introduce the families of complex networks that 
we consider in this work. Subsequently we define the no- 
tion of topological network entropy and apply it to study 
the structure of different complex network topologies. Fi- 
nally we briefly discuss relation with former works and 
state our conclusions. 

Network quantum states. Any graph Qa with n nodes 
is completely specified by its adjacency matrix A: a 2- 
dimensional array of size n, where each entry aij char- 
acterizes the connection between nodes i and j. The 
domain of aij determines the kind of graph one is con- 
sidering: directed, undirected, weighted or unweighted. 
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In this work we focus on undirected unweighted graphs, 
so called simple graphs, for which a^j £ {0, 1} (the 
same analysis can be naturally extended to directed and 
weighted graphs) . In particular if aij = it means there 
is no edge connecting the two nodes i and j, otherwise 
djj = 1. Sometimes it will be useful to refer also to non 
simple graphs with loops. Given the graph Ga {A ^ 0) 
we define the following bipartite quantum state 

n 

where \\A\\p '■= VTrA^A denotes the Frobenius norm of 
the matrix A and Hi = H-2 = C™. In the fixed local ba- 
sis {\i} : i — 1, . . . , n} we refer to \A) as a pure network 
state. It corresponds to the state of two n-level systems, 
or analogously to the state of two n-qubit systems where 
each subsystem of n qubits is constrained to the one- 
excitation manifold. The isomorphism class of a graph 
corresponds to the orbit of the permutation group on the 
adjacency matrix: PAP 1 , where P £ S n are permutation 
matrices. This implies that the adjacency matrix of iso- 
morphic graphs is unique up to permutations of rows and 
columns, and the same holds true for the bipartite states 
in EqUJ The reduced density matrix of the subsystem 
whose Hilbert space is Hi is given (in the given basis) by 

A /it 

PA =Tr 2 \A)(A\ = -^-eH 1 . (2) 

We refer to this reduced density matrix as the mixed 
network state. Notice that both definitions ((T|) and ([2]) 
do not rely on A being symmetric and therefore ex- 
tends immediately to oriented (and weighted) graphs. 
We are interested in the correlation properties between 
Hi and H2, quantified by the entanglement entropy 0]- 
Since this quantity depends only on the spectrum of the 
adjacency matrix, it is a property of the isomorphism 
class of the graph, i.e. isomorphic graphs will have the 
same entanglement entropy. Indeed, if P is an (orthog- 
onal) permutation matrix: \PAP f ) = P® 2 |A), where 
P = Yl7j=i ^».jl*)01- Namely isomorphic graphs give 
rise to locally equivalent network states. 

Before considering complex topologies, it is instructive 
to play with the simplest possible examples and try to 
characterize maximally and minimally entangled network 
states. Notice that the unnormalized bipartite state can 
be written as A<E> I\I), where \I) = J^i *) * s an unnor- 
malized maximally entangled state. The network corre- 
sponding to the state \I) consists of n nodes with loops, 
and by construction its adjacency matrix is the identity. 
Entanglement does not change under local unitary trans- 
formations Q, so in order to construct other maximally 
entangled network states we need to characterize all the 
adjacency matrices that correspond to unitary operators. 
It is easy to prove that the set of unitary adjacency ma- 
trices coincides with the set of permutation involutions, 



i.e. the permutation matrices that square to the identity. 
This is also consistent with the fact that the square of 
the adjacency matrix is the unnormalized totally mixed 
network state. Unitary adjacency matrices correspond 
to networks made of only loops or disconnected linked 
pairs of nodes. On the other hand factorized states (i.e. 
unentangled) correspond to complete graphs with loops. 

In the following we study properties of ensembles of 
random network states. The probability measure in the 
space of network states is the one induced by the mea- 
sure on the space of random networks, according to the 
construction in EqfT] and Eqf2] 

Complex networks. In order to make the paper self- 
contained, let us briefly introduce three network struc- 
tures that we will use in the following. The seminal paper 
of Erdos and Renyi in 1959 denned what is now the stan- 
dard example of a random network Q . The ER random 
graph model, denoted by Q^mi IS an ensemble of graphs 
where each element has n nodes and m edges. ER graphs 
are also related to so-called Gilbert models [||, denoted 
by Gn,pi where an edge between each pair of n nodes is 
present with a probability p. The Gilbert model is better 
suited for analytical investigations, while G®m graphs are 
numerically easier to study. In the thermodynamic limit, 
fixing the average degree q of a node, one can constrain 
the two models to be related by q = 2m /n — pn. If, 
for n —> oo, q/n — s- the network is said to be sparse. 
The other example of complex network that we consider 
is known as the Barabasi- Albert model [ltj, denoted by 
G BA , based on a growth process and a preferential at- 
tachment mechanism. The rationale is that nodes with 
higher degree acquire new nodes at higher rates than 
other lower-degree nodes. Nodes are added successively, 
and for each node a number d of edges are generated, 
with bias towards connections with higher-degree nodes. 
The distribution for the number of links emanating from 
a node is not Poissonian, like for ER graphs, but rather 
follows a power-law. Another way to model stochasticity 
in the connectivity pattern of a graph is by randomly de- 
stroying the periodicity of a regular lattice. This is the 
idea behind small- world networks, denoted by GpY > as 
proposed by Watts and Strogatz They can be cre- 
ated by randomly adding bonds to a regular one dimen- 
sional ring, this way building a superposition between 
regular lattices and classical random graphs. The prob- 
ability p according to which new bonds are added at 
random is a parameter characterizing the ensemble, and 
it allows to interpolate from regular graphs (p — 0) to 
ER random graphs (p — 1). The other parameter for 
this kind of networks is denoted with k, and it quantifies 
the number of next-nearest-neighbor links present in the 
original regular graph. For each of the above complex 
network ensembles we shall construct the associated en- 
semble of random network states denoted by ^ ER ^ BA 
and xjj sw , and we will consider scaling properties of the 
average entanglement entropies. 
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Topological network entanglement. We start by evalu- 
ating analytically the averaged Renyi entropy of network 
states in ip ER , the states associated to the ensemble of 
ER random graphs. The a-Renyi entropy of a state is 
defined by R a (p) = (1 — a) -1 log 2 Trp a . Using the defi- 
nition of pa given in Eq[2]we have 



R a (pa) 



log 2 TrA - a log 2 TrA 2 
1 — a 



(3) 



We are interested in the scaling in n of the average Renyi 
entropy R a . In order to provide an explicit expression for 
R a we use the fact that for each a there exists a constant 
C2a such that lim. n ^ > . 00 n~ 1 TrA 2a — c 2ct , i.e. for sparse 
ER graphs the thermodynamic limit of the moments of 
the graph spectrum exist and are finite [12j . Furthermore 
one can check numerically that the difference between 
the quenched average (1 — a) _1 log 2 p a and the annealed 
average (1 — a) _1 log 2/ o Q scales like n _1 . Putting this 
together we can write 



R a (pa) [n] = log 



,n + g(a) + 0(-). 

n 



(4) 



where g(a) := log 2 c 2 + (1 - a) 1 (log 2 c 2a - log 2 c 2 ) is 
a sub-leading O(l) term. This equation tells us that 
the Renyi entropy is almost maximal for any a. No- 
tice that, even though the logarithmic scaling for these 
network states is consistent with the one of general (Haar 



distributed) random states [14J, one could not predict a 



priori this behavior for the particular family of random 
states we introduced. Remarckably the sub-leading term 
contains information about the topology of the graphs. 
In fact, the term c 2a is directly related to the average 
number of closed paths of length 2a in the graph. In 
Fig we provide a numerical check of EqfJ] which sup- 
ports in particular the approximation of the quenched 
with the annealed average. The figure shows a perfect 
agreement between Eq|4] and the results of the simula- 
tion. It is tempting to extrapolate our analysis from 
the Renyi entropy to the von Neumann entanglement 
entropy, which is defined by S = lim a -n Ra- By con- 
struction it follows immediately the logarithmic scaling 
of the von Neumann entropy, while for the sub-leading 
term we have g e := lim Q ^i g(a) = log 2 c 2 - "f^ 2 * \ x =i 
and then, using the definition of c 2ct , one finds 



9e = log 2 C 2 



TrA 2 log 2 A 2 



TrA 2 



(5) 



It is interesting to observe that the second term in Eqj5] 
can be regarded as a sort of topological susceptibility of 
the given family of networks. In fact this term is equal 
c= i, and it tells us how the logarithm of the 



log£2i 

dx 



to 

rescaled averaged number of loops of length 2a changes 
as the length is changed continuously around a = 1. For 
this reason, and in view of its conceptual similarity with 
the topological entanglement entropy introduced in [l3[ 
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FIG. 2: (Color online) (a) The average Renyi entropy R a ob- 
tained numerically, and the analytical prediction (dashed and 
continuous lines) for ER quantum random states. The open 
square symbols corresponds to the a = 2 data points, while 
the open circles are for a — 3. (b) The average von Neumann 
entropy obtained numerically and the analytical prediction 
for the states corresponding to Gn,n- ( c ) <?e as a function of /3 
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(d) g e as a function of the rewiring probability 



p, for small world networks in Q^ 2 



we call the 0(1) quantity in Eq[5] topological network en- 
tanglement. In Fig(5}} it is shown a comparison between 
the analytical expression for the von Neumann entropy, 
obtained using Eq[5l and the empirical average of the en- 
tanglement entropy over different realizations. As can be 
seen from the figure there is a perfect overlap between the 
two. The ensemble Gnm ^ s parameterized by the num- 
ber of edges m. One can wonder about the behavior of 
g e as m varies in some particular interval, but in such a 
way that the graph is always sparse. Fig. shows the 
value of g e as a function of a parameter j3 G (0, 1] which 
is related to the number of edges by m = /3n log 2 n. The 
figure shows that for ER-graphs the greater the number 
of links the greater the entanglement. We now evaluate 
numerically g e for small worlds and preferential attach- 
ment networks. Considering first networks in we 
checked numerically that the scaling is logarithmic and 
that the prefactor is always 1. In Fig(5Ji we see the depen- 
dence of g e on p, the probability of rewiring edges. The 
figure shows that g e increases monotonically from regu- 
lar to more random graphs. This is consistent with the 
intuition that adding randomness to a graph increases 
its entropy, as measured by g e . From these results it 
is clear that the properties of the entanglement entropy 
provide information on the complex network structure, 
supporting its interpretation as a graph entropic quan- 
tity. Considering the ability to discriminate between dif- 
ferent network topologies, in Fig|3]we compare the scal- 
ing of the average von Neumann entropy of ip E R , i/jq 
and i\) BA random network states. The simulations show 
that, unlike for ER and small- world networks, the loga- 
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rithmic prefactor for ip BA states is slightly smaller than 
1. From the figure it is clear that the von Neumann 
entropy distinguishes different complex network ensem- 
bles. For sufficiently big networks the fluctuations due to 
disorder are strongly suppressed. On the one hand this 
is an indication of the robustnes of this graph-entropic 
measure; on the other hand it can also be useful from a 
computational point of view. In fact one has a very good 
estimate of the entanglement entropy already from few 
realizations. Hence, if the network is big enough the scal- 
ing analysis could in principle be done on one single real- 
ization, for example evaluating the entanglement entropy 
on sub-graphs of increasing size. Furthermore whenever 
computational efficiency is an issue we point out that 
instead of the Von Neumann entropy one could evaluate 
the so called single-copy entanglement (lim Q ,_ i . 00 R a ) 
for which efficient numerical techniques can be used 

Discussion and conclusions. In this letter we have ex- 
ploited a natural mapping from graphs to quantum bi- 
partite states and we have defined the entropy of a graph 
as the entanglement entropy of the associated quantum 
state. We have then used this quantum measure of corre- 
lations to study the structure of complex networks. The 
scaling of the entanglement entropy is logarithmic in the 
system dimension, and both the prefactor and the sub- 
leading 0(1) term (topological network entropy) can be 
used to characterize the network family and to distin- 
guish between different network topologies. In particular, 
we showed that the Barabasi- Albert model has a scaling 
behavior that differs significantly from the one of small- 
world and ER graphs. While these last two have a sim- 
ilar scaling, but still distinguishable comparing graphs 
with the same number of edges. This is consistent with 
the fact that small-world networks are mixture of regu- 
lar lattices and ER graphs. Furthermore, we provided an 
analytic expression, exact in the thermodynamic limit, 
for the averaged Renyi and von Neumann entropy asso- 
ciated to ER random graphs. It is desirable to achieve 
a clear and general understanding of the relations be- 
tween the quantum entropic measures we introduced and 
the standard graph-theoretic observables analyzed in the 
complex-network community. One would like also to gain 
a deeper insight into the measure concentration (large 
size convergence) properties of the various probabilistic 
objects we discussed for the different families of complex 
networks. While the primary goal of this paper has been 
to show how to use quantum tools to investigate complex 
networks it should be clear that also the converse task i.e., 
using properties of complex newtwork to study the novel 
class of random quantum states we introduced, is of in- 
terest on its own right. Moreover on the quantum side it 
is a challenge to find a consistent inverse mapping that 
allows one to associate to a general bipartite quantum 
state a specific network. Finally, one would like to devise 
efficient and physically feasible preparation schemes for 
the network quantum states we proposed. 
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FIG. 3: (Color online) The average von Neumann entropy 
obtained numerically for different quantum random states and 
the linear fit of the data points. 
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